Localized modes in nonlinear photonic kagome nanoribbons 
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We examine localization of light in nonlinear (Kerr) kagome lattices in the shape of narrow strips 
of varying width. For the narrowest ribbon, the band structure features a flat band leading to linear 
dynamical trapping of an initially localized excitation. We also find a geometry-induced bistability 
of the nonlinear modes as the width of the strip is changed. A crossover from one to two dimensions 
localization behavior is observed as the width is increased, attaining two-dimensional behavior for 
relatively narrow ribbons. 



The kagome lattice has persisted as an object of inter- 
est in several branches of physics since its appearance in 
1951 pQ. In condensed matter physics, kagome lattices are 
valued because they lead to spin frustration when the sys- 
tem contains antiferromagnetic interactions. Many prop- 
erties are still a matter of speculation, however, such as 
the nature of its magnetic ground state 0. In ultracold 
atom physics, several optical lattices have been used so 
far, most of which are primitive Bravais lattices. More 
recently, however, recent works have used non-standard 
optical lattices that contain a few-sites basis, such as the 
kagome lattice, hoping that the low energy degrees of 
freedom arising from the placing of a particle inside the 
unit cell, would bring richer dynamics and ordering[3]. 
Very recently, a magnetic field induced metal- insulator- 
transition has been predicted in kagome nanoribbons [2]. 
In photonics, kagome lattices have been recently pro- 
posed as a cladding in hollow-core polymer fibers, where 
it has been observed that the low overlap of the core 
modes and the microstructured material, plus the low 
density of states in the cladding, can lead to an improved 
core guidance [5]. The existence and stability of gap soli- 
tons and vortices in infinite, nonlinear kagome lattices 
have been examined in ref.[S]. In this work we focus on 
finite, ribbon-like kagome lattices and study the onset of 
localized modes and their stability as the width of the 
ribbon is changed, going from a quasi one-dimensional 
geometry to a quasi two-dimensional one. 

Let us consider a two-dimensional kagome photonic 
lattice in the form of a long ribbon, i.e., a N x M lattice 
of width w = M <C N (Fig.l). This type of waveguide 
array can be easily manufactured by means of the direct 
femtosecond laser- writing technology [7 . Another choice 
is the use of optical induction [S]. In the framework of the 
coupled-mode theory, the electric field E(r) propagating 
along the waveguides can be presented as a superposition 
of the waveguide modes, 

£(r)=££ n 0(r-n), (1) 

n 

where E n is the amplitude (in units of (Watt) 1 / 2 ) of the 
(single) guide mode <f>(r) centered on the site with lattice 



oooooooooooooo 
o o o o o o o 

oooooooooooooo 
oooooooo 

oooooooooooooo 
o o o o o o o 

oooooooooooooo 

Figure 1. Schematics of an optical kagome lattice nanoribbon 
of width w = 3. The ribbons are of indefinite extension in the 
horizontal direction, remaining finite in the vertical direction. 



number n = (711,712). The evolution equations for the 
modal amplitudes E n are: 

dE 

^+^Vn, m £ m + 7 |£ n | 2 £ n = 0, (2) 

m 

where n denotes the position of a guide center, z is the 
longitudinal distance (in meters), V n , m is the coupling 
between nearest-neighbor guides (in units of 1/mctcr) 
and 7 is the nonlinear coefficient (in units of l/(Watt 
x meter)), defined by 7 = lo^nij cA^, where ojq is the 
angular frequency of light, 712 is the nonlinear coefficient 
of the guide and A e Q is the effective area of the linear 
modes. The nonlinear parameter 7 is normalized to 1 for 
the focussing nonlinearity. 

Next, we analyze the stationary localized modes of 
Eq. (JlJ in the form E n (z) = E n exp(ij3z), where the am- 
plitudes E n satisfy the nonlinear difference equations, 

-PE a + J2 V n , m E m +j\E I1 \ 2 E Il ^0 (3) 

m 

For a given value of /3 the system is solved by a multi- 
dimensional Newton-Raphson method, using the anti- 
continuum (i.e., decoupled) limit as an initial condition. 
For the fundamental mode, we look for a localized so- 
lution with maximum at the center of the array and 
decaying along and across the ribbon. In order to vi- 
sualize the modes, we will use a generic gaussian mode 
<j){x,y) = exp[-(x 2 + y 2 )/a] in Eq.Q, with a = 0.15. 
as we will see, the results depend substantially upon the 
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ribbon width. Each mode family is characterized by its 
power content P = J2 n l-^«| 2 as a function of the propa- 
gation constant (i. 

The linear stability of each nonlinear mode is com- 
puted in the usual manner by defining a weakly per- 
turbed mode as E n (z) — (E n +u n (z)+iv n (z))exp(i/3z), 
where u n {z) and v n (z) are real. Their evolution equa- 
tions can be expressed in compact form by defining real 
vectors STJ = {«„} and 5V = {v n } and real matrices 
A = {A nm } = {<S„, m+ i + <5„,m-i + (~/3 + 37|£ , „| 2 )<5„, m 
andB = {B nm } = {6 ntm+1 +6 ntm -i + (-l3+j\E n \ 2 )6 n , m . 
The dynamical evolution of the perturbation can then be 
expressed as <5U + BA <5U = and SY + AB SY = 0, 
where an overdot denotes a derivative in z. Thus, the 
linear stability of the mode depends upon the eigenvalue 
spectra of AB and BA: If any of the real eigenvalues is 
negative, the mode is unstable; otherwise the solution is 
stable. 

Bulk modes: The narrowest ribbon has width w = 1 
and consists of a chain of hexagonal rings with a ba- 
sic cell containing 5 sites (Fig. 2). The linear (7 = 0) 
dispersion relation can be computed in closed form and 
consists of five bands: /3 = 2, j3 = ±y/2 + 2cos(fc), 
f3 = 1 ± 2cos(fc). The whole spectrum is gapless 

and one of the bands is completely flat, just as in the 
infinite lattice case. The existence of the flat band has 
implications for the diffusion of an initially localized ex- 
citation: If the initial profile has significant overlap with 
the modes on the flat band, then there will be a por- 
tion of the excitation that cannot propagate away from 
the initial site due to the vanishing group velocity of the 
band. As a consequence, there will be linear selftrapping 
at long evolution times. This localization effect is due to 
a destructive interference effect due to the inherent ge- 
ometry of the lattice. All this is borne out very nicely 
by a long-time computer simulation of the propagation 
of an initially localized pulse placed in the middle of the 
ribbon. Since the ribbon is really finite in the 'long' direc- 
tion, absorbing boundary conditions were used to avoid 
reflections from each end. The results are shown in Figj3] 
which shows the intensity of the field vs the evolution 
distance for the initial site and several other sites close 
and far from it. We see a localized optical profile around 
the initial guide, as predicted. 

For the nonlinear regime (7 7^ 0), Fig. 4 shows the 
power vs propagation constant curve for the ID funda- 
mental mode. It reaches all the way to the band edge, as 
is typical for ID lattices. 

We now increase the ribbon width to w — 2 and 
compute again the power vs propagation constant curve 
(Fig. 5). This time the curve features a stable and an un- 
stable regime, and does not reach the band edge. This is 
a two-dimensional-like behavior, which is surprising given 
the narrowness of the ribbon. For w = 3 (Fig. 6), we see 
still a different behavior: Now the ribbon behaves like 
a finite 2D lattice, characterized by a bistable power vs 
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Figure 2. Top: ID Kagome Ribbon. Bottom: Linear bands 
for a ID Kagome Ribbon. 
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Figure 3. Dynamical evolution of initially localized excitation. 
The plot shows the amplitude of the field as a function of 
longitudinal distance ('time'), at several sites away from the 
initial one. From top to bottom: At initial site, at a first 
nearest-neighbor, at a third nearest-neighbor and at an eighth 
nearest-neighbor . 



propagation constant curve. On each side of the bistable 
curve we have a stable mode whose width depends on 
its closeness to the band: It is wider (narrower) for the 
mode whose propagation constant is farther (closer) to 
the band edge. Separating these two stable modes, there 
is an unstable one. 

For w — 4 and higher (Fig. 7), we are back to the 
two-dimensional -like behavior, with a power vs prop- 
agation constant curve that is stable and decreases with 
decreasing propagation constant until its slope vanishes, 
followed by a unstable regime with negative slope, which 
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Figure 4. Left: Power vs propagation constant curve for 
the fundamental odd mode for a width w = 1 ribbon ("one- 
dimensional"). Vertical line marks the position of the band 
edge. Right: Intensity distribution for fundamental single- 
site mode marked with black circle. Shading is proportional 
to the intensity distribution of the mode power. 



is in agreement with the Vahitov-Kolokolov stability cri- 
terium for the fundamental mode[9]. 
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Figure 5. Left: Power vs propagation constant curve for the 
fundamental odd mode for a width w = 2 ribbon. Solid 
(dashed) portion denotes stable (unstable) regime. Vertical 
line marks the position of the band edge. Right: Intensity dis- 
tribution for fundamental single-site mode marked with black 
circle. Shading is proportional to the intensity distribution of 
the mode power. 

Surface modes: We compute nonlinear surface modes 
centered on a single waveguide at the very edge. A fi- 
nite ribbon sample features two different types of edge: 
Straight one along the 'long' direction and 'zig-zag' along 
the short direction (Fig.l). 

For a given ribbon width w, the power vs propagation 
constant curve for modes located on the straight edge 
exhibits the usual behavior already observed in other 
two-dimensional lattices, like the existence of a power 
threshold to generate the mode. This is observed for all 
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Figure 6. Left: Power vs propagation constant curve for the 
fundamental odd mode of width w = 3 ribbon. Solid (dashed) 
portion denotes stable (unstable) regime. Vertical line marks 
the position of the band edge. Right: Intensity distribution 
for fundamental single-site mode marked with black circle. 
Shading is proportional to the intensity distribution of the 
mode power. 
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Figure 7. Left: Power vs propagation constant curve for 
the fundamental odd mode of a width w = 4 ribbon. Solid 
(dashed) portion denotes stable (unstable) regime. Vertical 
line marks the position of the band edge. Right: Intensity dis- 
tribution for fundamental single-site mode marked with black 
circle. Shading is proportional to the intensity distribution of 
the mode power. 



w values examined, ranging from w — 1 up to w = 5. For 
widths greater than 4, the corresponding P vs (3 curves 
become indistinguishable from each other on the scale 
shown here. On the 'zig-zag' boundary there are three 
non-equivalent surface sites to consider, characterized by 
their different coordination numbers (2,3 and 4). Their 
power vs propagation constant curves are qualitatively 
similar. An example of such curves for w = 3 is shown 
in Fig. 8. The lowest threshold power in this case cor- 
responds to the mode centered on the zig-zag boundary 
with 2 nearest-neighbors. 

Thus, we have examined localization properties of 
two-dimensional photonic lattices consisting of kagome 
nanoribbons. For the quasi-one-dimensional case, and in 
the absence of nonlinearity, we found a flat band and 
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the ensuing linear localization of dynamical excitations, 
which agrees with the phenomenology found earlier for 
infinite kagome lattices. As the width of the ribbons 
increased, we found a crossover between one-and two- 
dimensional soliton behavior, characterized by the ap- 
pearance of a bistable regime, which constitutes an ex- 
ample of a geometry-induced bistability for solitons. This 
type of bistability could be generic of two-dimensional fi- 
nite lattices, as evidenced by similar results obtained for 
the graphene lattice [TP], The author is grateful to Y. S. 




°o° # o° 

o°o o°o 



5.5 - 



o o°o o 
o u o o o 



°o° 
o u o 



4.5 



4.5 



5.5 



°o° 
o u o 



4.5 



Figure 8. Left: Power vs propagation constant curve for sev- 
eral surface modes of a width w — 3 ribbon. Solid (dashed) 
portion denotes stable (unstable) regime. The insets show the 
vicinity of the array boundary and the the black circle marks 
the position of the mode center. 
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